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Abstract 

We pursue a number of analytical directions, motivated to some extent initially by the possibility 
of developing a methodology for formally proving or disproving a certain conjecture of quantum- 



theoretical relevance (quant-ph/0308037). It asserts that the 15-dimensional volume occupied by 



the separable two-qubit density matrices is {^/2 — 1) /3, as measured in terms of the statistical distin- 
guishability metric (four times the Bures or minimal monotone metric). Somewhat disappointingly, 
however, the several various analyses that we report, though we hope of independent/autonomous 
interest, appear to provide small indication of how to definitively resolve the conjecture. Among 
our studies here are ones of: (1) the Bures volumes of the two-dimensional sections of Bloch vectors 
for a number of the Jakobczyk-Siennicki two-qubit scenarios; (2) the structure of certain convex 
polytopes of separable density matrices; and (3) the diagonalization of 15 x 15 Bures metric tensors. 
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I. INTRODUCTION 

In a previous study we formulated — based on some combination of numerical and 
analytic evidence — a conjecture that the volume, as measured in terms of the statistical 
distinguishabihty (SD) metric j^, of the 15-dimensional convex set of 4 x 4 separable density 
matrices is 

= (V2 - l)/3 ^ 0.138071. (1) 

The numerator of ((T)), that is cr^^ = \/2 — 1 (alternatively, \/2 + 1 = \ — 1)), has been 
previously termed the "silver mean" P]. It is interesting to note that this SD conjecture 
involves the first three positive integers alone, while the "golden mean" jj] is defined as 
(V5 + l)/2 (or, alternatively as (^5 - l)/2 = 2/(v^+ 1)). 

A numerical integration based on two billion (separable and nonseparable) density ma- 
trices generated by a quasi-Monte Carlo (Faure-Tezuka) procedure had yielded an estimate 
of 0.137884 for '!]. (It now appears [3] that it is possible to significantly accelerate 
the MATHEMATICA program employed, so that we may soon be able to report results 
for substantially larger sample sizes.) In j^, an exact Bures probability of separability 
{^Bures/^Bure^"^^^^) GQual to aAg had becu obtained for both the q = I and q = \ two- 
qubit states (7|] inferred using the principle of maximum nonadditive [Tsallis] entropy 
and also for an additional /oty-dimensional scenario 

Since the Bures (minimal monotone j^) metric 
metric that is 

dBuresip, P + dpy = {I / 4:)dsD{p, P + dpf , (2) 

in the neighborhood of a density matrix p, the conjecture becomes equivalent to one 
that the Bures volume is 



6, sec. 2.2.11. 



KZes = 2-''(v^ - l)/3 ^ 4.2136 ■ 10"^ (3) 



KM 



(In [l|, we had also been led to a number of related conjectures, including one that Vj 
30V^^res' f*^^ two-qubit systems, where KM denotes the Kubo-Mori monotone metric.) 

In the present analysis, we undertake a line of research that hopefully — we, at least, 
initially thought — may contribute to formally proving/disproving this conjecture. This has 
seemed a quite formidable task, as a (naive "brute force") symbolic integration approach, 
along the lines of the high-dimensional numerical integration followed in jl|, appears to be 



far beyond present computing technology (cf. P]). In this previous work |])], we generated 
4x4 density matrices, which then had to be checked for possible separability. Here, our 



approach focuses at the outset on separah 
not directly enter the picture (cf. 
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e density matrices, and the nonseparable ones do 



3,0)- 

The convex domain (D) of 4 x 4 separable density matrices is 15-dimensional in nature. 
In sec. Ill Al we systematically generate 16 pure product 4x4 separable density matrices, 
convex combinations of these sixteen spanning some subset of D. (The entire domain D 
itself, however, is the closed convex hull of the set of all — not just 16 — product states 



|l4|.'l This approach can be interestingly contrasted (in sec. IIID|1 with that of Braunstein et 
al (l5[ , who used — in the 4x4 case — an overcomplete basis of 36 density matrices to "give 
a constructive proof that all mixed states of qubits in a sufficiently small neighborhood 
of the maximally mixed state are separable" . They allowed negative weights on their basis 
matrices in their analysis, thus enabling them to obtain entangled density matrices. We 
view the convex weights {wi,i = 1, . • . , 15) attached to the first of our 15 matrices as the 
parameters or variables of our problem. (Of course, then, we must have wiq = 1 — J^j^iWi 
bound or dependent parameter.) 
Proceeding onward, in sec. Ill Bl we investigate the possibility of using certain of the 
results in sec. Ill Al to obtain bounds on VBures- examine in sec. Ill CI the associated 
15 X 15 Bures metric tensor, and in sec. Ill Dl consider the use of an overcomplete basis. In 



sec. Ill E| we generate a polytope composed of 8 x 8 t/iree-qubit density matrices, while we 
obtain in sec. Ill FI the Bures metric tensor in diagonal form. We obtain and plot in sec. Ill (tI 



the Bures volume elements for t 



iree of the two-dimensional two-qubit scenarios presented 



by Jakobczyk and Siennicki jl6|, and derive the associated "Euclidean" probabilities of 
separability. 

II. ANALYSES 

A. Convex Polytope of 4 x 4 Density Matrices 

To generate the 16 extreme pure product basis matrices (cf. 12, Q]), we start with 
four 2x2 density matrices, the Bloch vectors of each of which extend to one of the four 
points ((f , f , f ), (— f , — f , f ), (— t;, f , — f ) and (f,— f,— f), where v = 1/VS) of an inscribed 
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tetrahedron (a Platonic solid). Then, we obtain 16 extreme points of by taking all possible 
■;ensor products of pairs of these four 2x2 density matrices. (We note that Schack and Caves 



19l sec. 2.4] mentioned the possibility of this construction, among others, observing that 
the "4 projectors for a tetrahedron are linearly independent, making the corresponding 
tetrahedral representation [of a 2 x 2 density matrix] unique".) 

These 16 extreme pure product density matrices have an interesting structure. Each 



one is a Bures distance of ^2 - 2/v^ ^ 0.919402 (and a Hilbert-Schmidt distance |20| of 
2/v^^ 1.1547) from six other matrices, and a Bures distance of 2/v^ ~ 1.1547 (and a HS 
distance of 4/3 ~ 1.3333) from the nine remaining matrices. (Using the trace of the matrix 
product of pairs of these matrices, as an alternative distance measure, we get the results 1/9 
[for six] and 1/3 [for nine].) The Bures distance between two density matrices {pi and P2), 
given by the formula. 



d 



Bures 



(pi,P2) = v/2-2Tr[(pfp2pf)^/^], (4) 



can only yield distances in the range [0, \/2]. The Hilbert-Schmidt distance is defined as 
eq. (2.3)] 

dHsipi,p2) = IIpi - P2\\hs = \/Tr[(pi - p2)^]. (5) 



Let us also observe that an explicit distance formula is presently availab 



"monotone" metric than the Bures, that is, the Wigner-Yanase metric 



21 



.e for just one other 
eq. (5.1)], namely. 



dwripi, p2f = 4arccos^ (tt[pI^^ pl^'^]^ . (6) 

(The class of all monotone metrics has, of course, a nondenumerable number of members.) 
In Fig. ^ we show a graph representation — in which only the links corresponding to the 



Bures distance \ 2- 2/V3 ^ 0.919402 are allowed — with the 16 extreme density matrices 

nn 

serving as the nodes. (The spectrum |22,|22l] of the associated (0,l)-adjacency matrix consists 
of the integers 6, -2 [nine-fold] and 2 [six-fold].) 

In Fig. 121 we show a graph representation — in which, now, only those edges corresponding 
to the Bures distance of 2/ -\/3 ~ 1.1547 are displayed. (Its spectrum consists of the integers 
9, 1 [nine-fold] and -3 [six-fold].) 
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FIG. 1: Graph representation of the 48 hnks between pairs of the 16 extreme pure product 



separable density matrices corresponding to a Bures distance equal to y2 — 2/-\/3 0.919402 
(and Hilbert-Schmidt distance of 2/\/3 « 1.1547). The 16 nodes are grouped into four quartets. 




FIG. 2: Graph representation of the 72 links between pairs of the 16 extreme pure product 
separable density matrices corresponding to a Bures distance equal to 2/\/3 ^ 1.1547 (and Hilbert- 
Schmidt distance of 4/3 ^ 1.33333.) 
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B. Volume calculations 



Additionally, the 16 extreme density matrices (and, it would seem any pure state, in 
general) are each a Bures distance of 1 (and a HS distance of \/3/2 ~ 0.866025) from the fully 
mixed state I4. (The WY-distance from each of the extreme 16 matrices to I4 is 27r/3.) Now, 
a 15-dimensional EuclideanhaJl of radius 1 has a volume equal to 2567r''/2027025 ~ 0.381443 
and one of radius ^3/2 ^ 0.866025, a volume equal to 27V37r73203200 ^ 0.044095. These 
Euclidean volumes are certainly larger than the Bures and HS volumes of the separable plus 
nonseparable 4x4 density matrices. These two latter values (precisely known only recently, 
making use of the theory of random matrices j3, eq. (4.12)] jsO, eq. (4.5)]) are, respectively. 



ysep+nonsep _ ^ nnnn574538- ysep+nonsep _ -.^^.^^^^-6 /„x 

^Bures " 165150720 ^ ^.0000574538, - 851350500 ^ ^'^^^^^ ' 

(Also, in ^], we had been led to conjecture that, in the two-qubit case, = 
64V^1J^"°"^'^^, where denotes the Kubo-Mori monotone metric.) This Bures volume, remark- 
ably, is the same as the 15-dimensional volume (" ') of a 15-dimensional 

n h n n 

Euclidean hemisphere (cf. j25[) with radius 1 [l|, |2J]. (In the numerical integration [l|, 
sec. 1.2] based on two billion sampled points, a very close estimate of 0.000574536 was 
obtained for V^^^esO Coupling this with the "silver mean" conjecture Q easily yields an 



implied conjecture that the Bures "probability of separability" 



;onieci 

" y 



IS 



psep _ ysZes _ (1680(^2 - 1) _ n^ooooQ (o\ 
^ Bures — ysep+nonsep — ^8 ~U.U/OOOOJ. ya) 

Bures 

The scalar curvature of the Bures metric at the fully mixed state J4 is known — as a specific 
case of more general formulas |27L Cor. 3] — to equal 570. We can, thus, attempt to apply 
the formula bs, p. 3676] (cf. pg, p. 936]), 
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for the volume of a geodesic ball centered at the n x n density matrix pq, where C(Po) is the 
scalar curvature at po- (The scalar curvature of the [monotone] Wigner-Yanase metric is a 
constant over the domain of n x n density matrices equal to |(n^ — l)(n^ — 2) [3^ Thm. 
7.2], which gives us "only" 52.5 for n = 4.) Taking the radius of the ball as r = 1, and 
C(po) = 570, n = 4 and Po = I, we get the result V{Bi{h)) = 5127r7883575 ^ -1.75015 



plus a higher-order term. (Note the negative sign! So, it appears that more exphcit terms in 
the asymptotic expansion @ are needed to try to avoid this. In his dissertation, A. Andai 
has such an expansion (with an O(r^) term), but considerably more differential-geometric 
information is needed beyond the scalar curvature to implement it (cf. |3,l|). 

If we fully mix, not the 16 extreme matrices, but only 15 of them in turn (the 15 nonzero 
Wj's all equalling 1/15), we find that the resulting 16 density matrices are all a Bures dis- 



tance of ^J2 - 2v/375 - 1/V5 ^ 0.0599422 (and a HS-distance of 1/(10^3) ^ 0.057735) 
from J4. A 15-dimensional Euclidean ball with this radius has a volume of only 
256(10 - V5- 2\/T5)^5/2^7(158361328125V5) ^ 1.76774 ■ 10-^^ (the comparable HS fig- 
ure being 1.00698 • 10^^^). Corresponding application of the formula © gives us a quite 
similar estimate of the Bures volume of 1.73225 ■ 10~^^ (plus a higher-order term) of our 
15-dimensional convex body of separable states. 

If we fully mix only 14 of the 16 (the 14 nonzero Wj's all equalling 1/14), then the resultant 
Bures distance from J4 is either 0.0810507 (for 72 of the possible pairs) or 0.0973228 (for 
the other 48 possible pairs). (The comparable HS-distances are -\/TT/42 ^ 0.0789673 and 



{^/5/3)/U 0.0922139.) 



Bures metric tensor 



Using our convex geometry/parameterization, we can further compute [32, Prop. 2] the 
15 X 15 Bures metric tensor (HsfijH) at I4. (This fully mixed state is obtained by setting 
Wi = 1/16 for all i, in forming the convex combination of the 16 extreme product basis 
matrices.) It has five eigenvalues equal to 17/s, where s = 98304 = 2^^ ■ 3, eight equal 



to 17/(3s) and a pair 17(31 ± v769)/(6s), having the approximate values 0.00169725 and 
0.000094224, being the roots of the quadratic equation 

(2^6 . 3^)x'^ - (2^1 ■ 3 ■ 7 ■ 31)x + 17^ = 0. (10) 

If we transform this tensor by the Jacobian (the determinant of which equals 2^^-\/2) cor- 
responding to returning from our choice of basis variables to that relying upon the naive 
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(11) 



parameterization of a two-qubit density matrix, say, 

/ an ai2 + 612^ ctis + hs^ + bui \ 

«12 — bui 0,22 0-23 + ^23^ «24 + &24^ 

ai3 - 6132 a23 - &23^ 033 a34 + 634? 

Vai4 - bui a24 - &24^ 034 - &34^ 1 - On - ^22 - ^33,/ 

we can obtain the invariant trace of the Bures metric tensor (cf. js^), which is 255/2^^. 
The square root of the determinant of the transformed tensor gives us the invariant volume 
element at h- This we computed to equal 2-120177^17/2 ^ 9.00021 ■ lO^^s. 

If we allow the weight (wi) assigned to the first, say, of the 16 matrices to vary between 
and 1/8, and the weight {wiq = 1/8 — Wi) compensatingly likewise, keeping the fourteen 
other weights as 1/16, the determinant of the 15 x 15 Bures metric tensor is a (high) 80- 
degree polynomial in wi, but the graph (Fig. of its square root multiplied by the Jacobian 
determinant (that is, the corresponding volume element l^fjjl^/^) seems quite well-behaved — 
clearly peaking at wi = wiq = 1/16, that is at J4. The integral of this function over [0,1/8] 
is 9.38545 ■ lO^^^ — a result certainly not inconsistent with our "silver mean" conjecture Q 
for the total Bures volume V^^^g^ the full 15-dimensional convex set of 4 x 4 separable 

states. The Bures distance between the two extremes on this plot (that is, one state with 
Wi = 0, 1^16 = 1/8 and one with wi = 1/8, wiq = 0, while all other w's remain fixed at 1/16) 
is (6 — -\/34)/6 ~ 0.0281747, considerably larger than the indicated integral between these 
two endpoints. In Fig. ^ we show the same form of plot as Fig. El but now using the trace 
of the 15 X 15 metric tensor rather than the volume element. 

We were not able to compute determinants of the 15 x 15 metric tensor when more than 
one of the convex weights was allowed to freely vary. (We were also unable to express a Bell 
state as a linear combination — of course, requiring negative weights — of the 16 extreme 
pure product density matrices.) 



D. Overcomplete Basis 

If we were to use in the two-qubit case, not the 16 extreme density matrices employed 
above, but the 36 extreme density matrices (forming an overcomplete basis), applied by 
Braunstein et al then we find that of the 35-36/2 = 630 possible distinct unordered 
pairs of density matrices, 144 correspond to a Bures distance of ^2-^2, 198 to V2 and 
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FIG. 3: Plot of volume element of 15 x 15 Bures metric tensor in the neighborhood of the fully 
mixed state /4, as a function of wi over the range [0,1/8], with Wi = 1/16 for i = 2, . . . , 15 and 
wiQ = 1/8 — wi. For wi = wiQ = 1/16, one gets I4. The function peaks at 14, having the exact 
value 2-120177^x7/2 PS 9.00021 • lO'^^ there. 
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FIG. 4: Plot of trace of 15 x 15 Bures metric tensor in the neighborhood of the fully mixed state 
/4, as a function of wi over the range [0,1/8], with Wi = 1/16 for i = 2, . . . , 15 and wiq = 1/8 — Wi. 
For wi = wie = 1/16, one gets I4. The function peaks at wi = 1/16 with the value 255/2^^ ~ 
0.00389099 there. 

288 to 1. (Thus, as opposed to the analysis based on the 16 matrices, there are three, not 
two different sets of internodal distances.) In Fig. we represent those 144 links with a 
Bures distance of 
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FIG. 5: Graph representation of the 144 hnks between pairs of the 36 extreme pure product 
separable density matrices — based on the overcomplete "Pauh" basis of Braunstein et al for 
the two-qubit case — corresponding to a Bures distance equal to \/ 2 — \f2. The nodes are grouped 
into sextets. 

E. Convex Polytope of 8 x 8 Density Matrices 

We also extended the general approach above to the t/iree-qubit case, by taking tensor 
products of all possible triplets of the same four 2x2 density matrices with their Bloch 
vectors reaching the points of the inscribed tetrahedron. (The extension to the quhit- qutrit 
case seems more problematical, since there seems to be no natural immediate analogue of 
the tetrahedron to employ for choosing nine 3x3 density matrices to tensor product with 
the four 2x2 density matrices we have been utilizing.) This gives us 64 extreme product 
triseparable states. There are, then, 63 ■ 64/2 = 2016 distinct pairs of intermatrix Bures 
distances. Of these, 288 are y^2 - 2/VS ^ 0.919402, 864 are 2/^3 ^ 1.1547 and another 
864 are ^2 - 2/(3^3) ^ 1.27087. (The corresponding HS-distances are 2/VS ^ 1.1547, 
4/3 ^ 1.33333 and 2^/13/3/3 ^ 1.38778.) 

In Fig. ini we show a graph representation — in which only the edges corresponding to 
the Bures distance ^2 — 2/-\/3 ~ 0.919402 are allowed — with the 64 extreme density 
matrices serving as nodes. (The spectrum of the associated (0,l)-adjacency matrix consists 
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FIG. 6: Graph representation of the 288 hnks between pairs of the 64 extreme pure product 
triseparable density matrices having Bures distance equal to \J'^ — 2/\/3 ~ 0.919402 (and Hilbert- 
Schmidt distance equal to 2/^/?, k, 1.1547). The 64 nodes fall into 4 groups of 16. 

of the integers 9, 5 [nine- fold], -3 [27-fold] and 1 [27- fold]. For the graph based on the 
Bures distance 2/-\/3, the spectrum has 27, -5 [27- fold], and 3 [36-fold] and that based on 
y2^^2/(3Vi), 27, -9 [nine-fold], 3 [27-fold] and -1 [27-fold].) Graph representations with 
the other two possible links are perhaps too dense visually for a meaningful representation. 

The 64 extreme matrices are all a Bures distance of ^2 — 1/ \/2 ^ 1.13705 (and a HS- 
distance of ^7/2/2 ^ 0.935414) from the fully mixed state 1%. Euclidean 3 5- dimensional 
balls of these radii have volumes of 0.0000299764 and 3.23277 ■ 10~^, respectively. (The 
formula Q — taking into account that the scalar curvature at /§ equals 9954 27] — gives 
us a "volume" 1.03199 ■ 10"^'*, but with a minus sign.) The trace of the associated 63 x 63 
Bures metric tensor was found at Is to equal 7/2^^. The square root of the determinant 
of the tensor equals (2^^^^3''^\/2)"^. (Multiplication by the determinant of the Jacobian of 
the transformation back to the naive parameterization of the 8x8 density matrices (cf. 
(II 1|) ) would give us the volume element at /§.) The spectrum of the tensor consists of 26 
eigenvalues equal to l/(2'^^3^), 26 equal to l/(2'^^3^) and 8 equal to (l/(2'^^3) — plus three 
(1.90525-10-22, 1.3992-10-24 and 6.10756-10-^4) — that are the roots of the cubic equation, 
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(2219 . 36)^3 _ (2145 . 33 . loi)^^ + (2^9 ■ 3 ■ 499)a; -1 = 0. (12) 
We have not, however, yet been able to compute the associated eigenvectors. 

F. Diagonal 15 x 15 Bures Metric Tensor 

We were not only able — using MATHEMATICA — to find the (four distinct) eigen- 
values, given above in sec. Ill C\ using the parameterization based on the 16 extreme pure 
product density matrices, of the 15 x 15 Bures metric tensor at J4, but also the associated 
set of orthonormal eigenvectors. Using the components of these eigenvectors as weights, we 
performed a linear transformation from the variables Wi to a new set, call them Vi, so that 



the metric tensor was now diagonal (cf. '34l. l35l|). Additionally, at J4, only two of the Vi^s 



were nonzero — let us denote them vi and f 7 (as this is how they arise in the analysis) — 
and these were precisely the ones associated with the paired eigenvalues that are the roots 
of the quadratic equation (|TIHl . So, at I4, all fj's are zero, but for two that are the square 
roots of the roots of the quadratic equation 

(2^ ■ 769)x^ - (2 ■ 3 ■ 5 ■ 769)x - 3^ = 0. (13) 



These square roots are y ^ ± 512^^759 ' '^^^i^S the approximate values of 0.239581 and 
0.0345646. (The larger value is associated with the eigenvector for the larger root of (fTUI) .) 
The 15 4 X 4 density matrices weighted by the 15 Wj's in the new expansion of a 4 x 4 
density matrix were all traceless 4x4 Hermitian matrices. (They were mutually orthogonal, 
but their self-inner products did not all have the same value, though certainly they could 
be rescaled so that they do. If we do, in fact, rescale so that all self- inner products equal 



2, then we have that fi,f7 = 3 ± -^^O There was also a "constant" term (weighted by 
no parameter Vi) in this expansion, which was the density matrix of a pure state. We can 
write this density matrix as \j){j\, where 

\J) = , ' hi + ^)(l + v^), hi + + v^), !)• (14) 

^6(2 + ^3) 1 + ^3 2 2 2 2 



Just as we completed the analysis immediately above, we noticed the appearance of a 
preprint of Kimura and Kossakowski (cf. j^^) in which the (somewhat similar) density 
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operator expansion 

P = ^lN + l^Sr\TrpX.^X, (15) 

is proposed. Here the A's {i = 1, . . . , A^^ — 1) are orthogonal generators of SU{N) which 
satisfy 

(i)Al = A„ (ii)Tr[A,] = 0, (iii)Tr[A,A,] = 26,,. (16) 

We, then, reran our computer programs using this parameterization (employing the standard 
generators of SU{A)). We found that for the fully mixed state I4 (where all (Aj) = Tr[pAj] 
are equal to zero), not only that the 15 x 15 Bures metric tensor was diagonal, but also 
that all its diagonal entries (eigenvalues) were equal to a common value, that is, 17/262144, 
where 262144 = 2^*^. The square root of the determinant of the tensor was, then, equal to 
(17/262144)15/^ which upon multiplication by the determinant 2^'^\/2 of the Jacobian for 
the change-of- variables back to the naive parameterization (lllj) of the 4x4 density matrices, 
gives us again the invariant volume element 2~i^°17^a/17/2 9.00021 ■ 10"^*^. 

Also, when we let one parameter vary freely ((Ai), say) from 0, the Bures metric tensor 
remained diagonal, unlike the situation in our earlier analysis. However, it appeared that if 
we were to let all 15 of the parameters vary freely, then the 15 x 15 Bures metric tensor is not 
diagonal. (We generated the [very large] symbolic entries of this tensor and then substituted 
random values of the 15 parameters, in reaching this conclusion.) This has the practical, as 
well as theoretical significance that one can readily compute the determinant of a symbolic 
multivariate diagonal 15 x 15 matrix, but certainly not one with many off-diagonal nonzero 
entries. 

The 15 X 15 Bures metric tensor obtained using the naive parametrization (jlip is itself 
quasi- diagonal, with the only off-diagonal structure being the 3x3 submatrix pertaining to 
the diagonal entries 011,022,033. The six off-diagonal nonzero entries are 2"^^ and the 15 
diagonal entries are 2^^^. The volume element (square root of the determinant) at J4 is, as 
mentioned previously, 2^"'^^°17'^a/17/2. 



G. Jakobczyk-Siennicki Two-Qubit Analyses 

Jakobczyk and Siennicki studied — using the same SU{4) parameterization (jl5p as 
employed by Kimura and Kossakowski 3^ (which we used to obtain a diagonal equal- entry 
15 X 15 Bures metric tensor [sec. IIIFj ) — two-dimensional sections of a set of Bloch vectors 
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FIG. 7: Bures volume element over the triangular domain of two-parameter 4x4 (separable and 
nonseparable) density matrices, denoted C by Jakobczyk and Sennini Fig. 1]. The Bures 
volume was computed to be 6.43885 • 10^'^. The fully mixed state I4 is at the origin. 

corresponding to two-qubit density matrices. In an interesting set of diagrams, they depicted 



both the domains of these 4x4 density matrices 



16l Fig. 1] and their convex subsets of 



separable states [16|, Fig. 2] (in those specific cases in which there were nontrivial inclusions). 

We have, first, analyzed the case they denote by the letter C, in which only the sixth and 
fifteenth SU{4) generators are assigned nonzero weight in the density matrix expansion (|15|). 
We found that the total Bures volume (of the indicated triangular domain) was 3.27995-10"^. 
Then, using further results of theirs, we computed the Bures volume of the separable subset 
to be 6.43885 ■ 10"''. Thus, the probabihty of separability (the ratio of these two quantities) 
j2^,Q for this particular scenario (labeled CK') was 0.196309. (Viewing the corresponding 
illustration of the separable and nonseparable regions in |l6l . Fig. 2] simply naively as a 
Euclidean diagram, we obtain a "Euclidean probability of separability" that is much higher, 
being the ratio of 1/^6 + v^7r/9 to 2^/73^/^ that is, (9 + 27rV3)/24 ^ 0.82845. Also, we 
were not able to compute enough digits of accuracy to meaningfully address the question 
of whether these quantities might possibly correspond to analytically exact expressions.) In 
Fig. 13 we show the Bures volume element over this two-dimensional triangular domain of 
4x4 density matrices. (The fully mixed state J4 corresponds to the origin.) 

In Fig.lHlwe analogously show the Bures volume element for the two-dimensional scenario 
Jakobczyk and Siennicki denote by the letter G (in which only the sixth and eighth SU{4) 
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0.5 



FIG. 8: Bures volume element over the elliptical domain of two-parameter 4x4 (separable and 
nonseparable) density matrices, denoted G by Jakobczyk and Sennini Fig. 1]. 

generators are allowed to bear nonzero weight). As a Euclidean figure, the area of this 
elliptical domain is (97ry372)/32 ^ 1.08215, while the area of the separable subdomain is 
(9(2 + 7r)v/372)/64 ^ 0.885535. The naive Euclidean probability of separability is, then, 
simply (2 + 7r)/(27r) ^ 0.81831. Our attempts to compute with reasonable confidence the 
corresponding Bures volumes (and, thus, the Bures probability of separability) were impeded 
by considerable MATHEMATICA numerical integration difficulties we have not yet resolved. 

In Fig. El we further show the Bures volume element for the two-dimensional scenario 
Jakobczyk and Siennicki denote by the letter E (in which only the third and nineth SU{A) 
generators are allowed to bear nonzero weight). As a Euclidean figure, the area of this 
hyperbolic domain is 2v^/3 ~ 0.942809, while the area of the separable subdomain is simply 
2/3 ~ 0.666667. The naive Euclidean probability of separability is, then, 1/v^ ^ 0.707107. 
Our attempts to compute with reasonable confidence the corresponding Bures volumes (and, 
thus, the Bures probability of separability) were again impeded by numerical integration 
problems. 



III. CONCLUDING REMARKS 



Let us point out a number of papers 



40 



4l| (see also j43. Chap. 6]), in which 



the volumes of non-Euclidean regular polytopes — such as we have encountered in this 
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FIG. 9: Bures volume element over the hyperbolic domain of two-parameter 4x4 (separable and 
nonseparable) density matrices, denoted E by Jakobczyk and Sennini 16, Fig. 1]. 



study (sees. Ill AIIII E|) — are computed. For an (Euler ang 
density matrices different from those discussed above, see 



e) parameterization of the 4x4 
J43|. (This parameterization was 
used in the numerical integration in which the "silver mean" conjectures ((0), ©) were 
formulated, as the domain of possible values of the parameters is readily expressible as a 
15-dimensional hyperrectangle.) 

Let us conclude on something of a discouraging note, however, in so far as we had been 
hoping initially to be able to report some progress or conceptual breakthrough in formally re- 
solving our "silver mean conjectures" (((T)), Q), as to the statistical distinguishability/Bures 
volumes of the 15-dimensional convex set of separable two-qubit states. (We continue to 
pu..ue the po^Mity ot developing companion coniectu.e= fo. the 35-d.men=io„al convex =et 
of quhit- qutrit states '5'], in terms of the Bures and a number of other monotone metrics, as 
well as the Riemannian, but non-monotone Hilbert-Schmidt metric M\-) The analyses 
presented above — though perhaps of interest from a number of perspectives — appear to 
help little in this regard. 



One possible further approach we have not yet investigated in any depth is the application 



43], seeing that the Bures metric does, in fact, serve 



of the concept of mmzma/ volume 

as the minimal monotone metric S|- (For a smooth manifold M, the minimal volume of M 
is defined to be the greatest lower bound of the total volumes of M with respect to complete 
Riemannian metrics, the sectional curvatures of which are bounded above in absolute value 



16 



by 1 j48|. We note that /ower bounds on the scalar curvatuie are available from the work of 



Dittmann 



H ) In 3 



|, demonstrating a conjecture of Gromov, the minimal volume of R 



2 



the infinite Euclidean plane) was shown to be (An exposition of this result is given in 
.) In |42| the value of was obtained for a certain supremum of volumes. 
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